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O ' Abstract 

a^ ■ 

0^ ■ We classify (l+3)-dimensional Fokker-Planck equations with a 

i constant diagonal diffusion matrix that are solvable by the method 

Oh| of separation of variables. As a result, we get possible forms of the 

■ drift coefficients Bi (x) , B2 (x) , B-^ (x) providing separability of the cor- 

, responding Fokker-Planck equations and carry out variable separation 

in the latter. It is established, in particular, that the necessary condi- 
tion for the Fokker-Planck equation to be separable is that the drift 
coefficients B{x) must be linear. We also find the necessary condition 
r> I for R-separability of the Fokker-Planck equation. Furthermore, exact 

' solutions of the Fokker-Planck equation with separated variables are 

constructed. 



I. Introduction 

The diffusion processes play an important role in different fields of physics, 
chemistry and biology |1|, |^ and have very broad applications in technics. For 
example, they play a decisive role in electronics That is why, analytical 
solutions of diffusion equations were, are and will be of great importance for 
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applications, since they provide new insights into the nature of the diffusion 
processes described by the equations in question. As an example, lets us men- 
tion the famous Black-Scholes models, whose success is based, in particular, 
on an analytical solution constructed in an explicit form. 

Unfortunately, the diffusion (Fokker-Planck) equations, that are used in 
applications, have variable coefficients and cannot be integrated by the stan- 
dard Fourier transform method. In fact, the only available efficient way for 
constructing analytical solutions of partial differential equations (PDEs) with 
variable coefficients is the method of separation of variables. 

The simplest and most widely used in applications is the case of constant 
diffusion matrix. Therefore the principal object of the study in the present 
paper is a problem of separation of variables in the Fokker-Planck equation 
(FPE) with a constant diagonal diffusion matrix 

Ut + Au + {Ba{x)u),^ = 0, (1) 

where B{x) = {Bi{x) , B2{x) , B^{x)) is the drift velocity vector. Here u = 
u{t,x) and Bi{x), i = 1,2,3 are smooth real-valued functions. Hereafter, 
the subscript Xa implies partial differentiation and, moreover, summation 
over the repeated Latin indices from 1 to 3 is understood. The matrix of the 
constant diffusion coefficients is reduced to the unit matrix by proper simple 
transformation. 

In the present paper we solve the problem of variable separation in FPE 
(|l]) into second-order ordinary differential equations in a sense that we obtain 
possible forms of the drift coefficients Bi{x),B2{x),B^{x) providing separa- 
bility of (|1]). Furthermore, we construct inequivalent coordinate systems en- 
abling to separate variables in the corresponding FPEs and carry out variable 
separation. 

The separability criteria for the one-dimensional FPE have been obtained 
in §. 

The problem of variable separation in the three-dimensional FPE was 
considered in a restricted sense by Sukhomlin in |]^. He used the symmetry 
approach which is based on the well-known fact that a solution with sep- 
arated variables is a common eigenfunction of three first- or second-order 
differential operators which commute with each other and with the opera- 
tor of the equation under consideration (see, and the references therein). 
Sukhomlin has obtained some drift coefficients providing separability of FPE 
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(|I]) and carried out separation of variables in the latter. But his results are 
far from being complete and systematic. 

Our analysis is based on the direct approach to variable separation in 
linear PDEs suggested in 0-0. It has been successfully applied to solving 
variable separation problem in the wave and Schrodinger equations [p|-|ll2 
with variable coefficients. 



II. Separation of variables in the Fokker— Planck 
equation 

Let us formulate briefly the algorithm of variable separation in FPE (|^) 
following 0. 

We say that FPE (|I|) is separable in a coordinate system t, uja = uja{t, x), 
a = 1, 2, 3 if the separation Ansatz 

3 

U{t, X) = ipo{t) Yl ipa {uJa{t, x) , X) (2) 
a=l 

reduces PDE (|l]) to four ordinary differential equations for the functions 
(ff,, in = 0, 1,2,3) 

V^O = Uo{t, ipo; A), ip" = Ua{uJa, (fa, <~p'a\ X). (3) 

Here f/o, . . . , f/3 are some smooth functions of the indicated variables, A = 
(Ai, A2, A3) G A = {an open domain in R^} are separation constants (spectral 
parameters, eigenvalues) and, what is more. 



rank 



3 3 



= 3. (4) 

/i=0 a=l 



The above condition secures essential dependence of a solution with separated 
variables on the separation constants A. 

The principal steps of the procedure of variable separation in FPE (|^) 
are as follows 

1. We insert the Ansatz (0) into FPE and express the derivatives (^q, 

5 V's in terms of functions y^o, v^i, (p^-, V'3, V^'i, ^P^i V^s using equations 
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2. We regard y^o, fi, ip2, ^3, V^'i, '^'2, V^s, ^i, ^2, A3 as the new indepen- 
dent variables yi,. . ., yiQ. As the functions ui, 002, are independent 
of the variables yi, . . . ,yio, we can split by these variables and get an 
over-determined system of nonlinear partial differential equations for 
unknown functions Ui, UJ2, UJ3. 

3. After solving the above system we get an exhaustive description of 
coordinate systems providing separability of FPE. 

Having performed the first two steps of the above algorithm we arrive at 
the conclusion that the separation equations (^) are linear both in y^o, • • • , V's 
and Ai, A2, A3. 

Next, we introduce an equivalence relation S on the set of all coordinate 
systems providing separability of FPE. We say that two coordinate systems 
t,ui,uj2,uj3 and t,u;i, 002,^^3 are equivalent if the corresponding Ansatzes (|^) 
are transformed one into another by the invertible transformations of the 
form 

t^i=fo{t), uji = fi{uji), (5) 

where /o, . . . , /s are some smooth functions and i = 1,2,3. These equivalent 
coordinate systems give rise to the same solution with separated variables, 
therefore we shall not distinguish between them. The equivalence relation (|^) 
splits the set of all possible coordinate systems into equivalence classes. In a 
sequel, when presenting the lists of coordinate systems enabling us to separate 
variables in FPE we will give only one representative for each equivalence 
class. 

Following [|1^ we choose the reduced equations (§) to be 

(P'q = (To(t) - Ti(t)Ai) V^O, = {Fao{uJa) + Fai{uJa)Xi) '-Pa, (6) 

where Tq, Tj, Fao, Fai are some smooth functions of the indicated variables, 
a = 1,2, 3. With this remark the system of nonlinear PDEs for unknown 
functions tUi, uj2., UJ3 takes the form 

1,2,3; (7) 
a = 1,2, 3; (8) 



duoi dujj 



= 0, i,i 



Tait), 
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B,^ + ^ + Aua = 0, a = l,2,3; (9) 

Thus the problem of variable separation in FPE reduces to integrating 
system of ten nonlinear PDEs for three functions. What is more, some coeffi- 
cients are arbitrary functions, which should be determined while integrating 
equations ([7p-(^). We have succeeded in constructing their general solution 
which yields, in particular, all possible functions Bi{x), B2{x), B^i^x) such 
that FPE (|l]) is solvable by the method of separation of variables. 

The system of equations (|^, (H) has been integrated in [Q. 

Lemma 1 The general solution uo = uj{t, x) of system of partial differential 
equations ^ is given implicitly by the following formulae: 

x = T{t)H{t)z{uj) + w{t). (11) 

Here T{t) is the time- dependent 3x3 orthogonal matrix: 

cos a cos /3 — sin a sin (3 cos 7 
T{t) = I sin a cos /? + cos a sin /5 cos 7 — ^ (12) 

sin P sin 7 

— cos a sin /5 — sin a cos /3 cos 7 sin a sin 7 

— sin a sin /? + cos a cos /? cos 7 — cos a sin 7 
cos/? sin 7 cos 7 



a,/3,7 being arbitrary smooth functions oft; z = z{u) is given by one of the 
eleven formulae 

1. Cartesian coordinate system 

Zl = UJi, Z2 = UJ2, z^ = u^. 

2. Cylindrical coordinate system 

Zi = e*^^ cosci;25 ^2 = e'^^ sin 1^25 2:3 = lj^. 

3. Parabolic cylindrical coordinate system 

Zi = {ujI-UjI)/2, Z2 = UJlUJ2, ^3=^3- 

4. Elliptic cylindrical coordinate system 
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Zi = a cosh uji cos 002, Z2 = a sinh ooi sin 002, Z3 = 003. 

5. Spherical coordinate system 

Zi = u;f ^ sech U2 cos u^, Z2 = uj^^ sech U2 sin tus, 
= cuf ^ tanhti;2. 

6. Prolate spheroidal coordinate system 

zi = a csch ui sech u;2 cos cus, 2:2 = a csch cui sech U2 sin ^3, 

2:3 = acothtui tanhc<j2. (13) 

7. Oblate spheroidal coordinate system 

Zi = a sec ciJi sech uj2 cos ^^3, ^2 = a sec loi sech a;2 sin a;3, 
Z3 = atancui tanhci;2. 

8. Parabolic coordinate system 

= e'^i+'^2 C0SCJ3, Z2 = e'^i+'^2 sinu;3, 
= (e^-i -e'"^)/2. 

9. Paraboloidal coordinate system 

zi = 2a cosh LJi cos ^^2 sinh ^3, Z2 = 2a sinh cui sin U2 cosh u;3, 
^3 = a(cosh2ci;i + cos2ci;2 — cosh2ti;3)/2. 

10. Ellipsoidal coordinate system 

Zi = ik~^ {k')~^dii{iUi, k) dn(co'2, k) dn(co'3, k), 
Z2 = —k{k')^^cn{uji, fc) cn(co'2, ^) cn(co'3, k), 
Zs = ksn{ui, k)sn{u2, k)sn{u!s, k). 

11. Conical coordinate system 

Zi = uj^^{k/)^'^dn{u2, fc) dn(u;3, k), 
Z2 = iuji^k(k')~^cn(uj2, A;) cn(u;3, k), 
Z3 = Ui^ksn{tU2, k) sn(ti;3, k); 
H(t) is the 3x3 diagonal matrix 

hit) 

Hit) = I h2it) I , (14) 








where 



(a) hiit), h2it), h2it) are arbitrary smooth functions for the completely split 
coordinate system (case 1 from jjiTSj)), 
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(b) hi{t) = h2{t), hi{t), hslt) being arbitrary smooth functions, for the par- 
tially split coordinate systems (cases 2-4 from ^TBj)), 

(c) hi(t) = h2{t) = h^{t), hi{t) being an arbitrary smooth function, for 
non-split coordinate systems (cases 5-11 from (fZ^jj 



and w{t) stands for the vector-column whose entries Wi{t),W2{t),W3{t) are 
arbitrary smooth functions oft. 

Here we use the usual notations for the trigonometric, hyperbohc and Jacobi 
elhptic functions, k {0 < k < 1) being the modulus of the latter and k' = 
(1 — A;^)^/^. To obtain real values for Zi,Z2, for the ellipsoidal coordinates 
(system 10) we choose ui real, uj2 complex such that Re uj2 = K, and 
complex such that Im uj^ = K', where K, K' are defined by 



7r/2 



K{k) 



To cover all real values of 2:1, 2:2, 2:3 once, it is sufficient to let uji vary in the 
interval [—-ft', K], 002 vary in [K — iK\ K + zi^'] (parallel to imaginary axis), 
and tus vary in [—K-\-iK\K + iK']. For the conical coordinates (system 11) 
uji,uj2,ui3 have the range < tui, —2K < uj2 < 2K, K < < K -\- 2iK'. 



For more details about elliptic functions, see |14 



Moreover we have obtained the explicit forms of the functions F^, [i,j = 
1,2,3) for each class of functions z = z{uj) given in ([T3D . The results are 
presented below in the form of 3 x 3 Stackel matrices |]T3| J-'i, . . . , J-'u, whose 
(i, j)th entry is the corresponding function Fij{uji). 



^1 




a cosh uoi 1 



a? cos^ UJ2 








-^1 

,-2 











cosh 002 — 1 



1 



7 



^6 



( sinh ^ liJ\ — sinh ^cji — 1 

COsh"*^ LO2 COsh~^ ijJ2 — 1 



v 







a} cos ^ uj\ 







- cos CJi 

— cosh~^ UJ2 cosh^^ uj^ 




(15) 



g4a;i _p2u)x 



e 

V 



4a;2 



-e 




^1 



10 



-^11 



cosh 2ti;i —a cosh 2uji 
— cos^ 2^72 a cos 2ci;2 
cosh 2uj'\ a cosh 2ct;cj 



/ sn^(u;i,A;) sn^(u;i,A;) 



cjf^ -c^r^ 

-khYi\uj2,k) 1 
-A;W(tU3,A;) 1 




We have also got the expressions for Ti(t), T2(t), T3(t) in terms of /ii(t), 
/i2(t), /i3(t): 



1. = hT\ 

2-4. Ti = /i^^ 
5 - 11. Ti = /i^^ 



z = 1,2,3; 
T2 = 0, T3 : 
T2 = T3 = 0. 



(16) 



In view of the above it is not difficult to integrate the remaining equations 
from the system under study. 

Note that we have chosen the coordinate systems uji,uj2,uj3 with the use 
of the equivalence relation S (H) in such a way that the relations 



AuJa = 0, 



1,2,3 



(17) 
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hold for all the cases 1-11 in (]T3|). Solving (^) with respect to Bj{x),i = 1, 2, 3 
we get (see, also [Q) 

B{x) = M{t){x-w)+'iv. (18) 
Here we use the designation 

M{t) = f{t)T-\t) + T{t)H{t)H-\t)T-\t), (19) 



where T(t), H{t) are variable 3x3 matrices defined by formulae ([T^) and 
([liD , correspondingly, w = (i(;i(t), ti?2(t), ti'3(t))"^ and the dot over a symbol 
means differentiation with respect to t. 

As the functions Bi, B2, B^ are independent of t, it follows from (|T8|) that 

-B(f) = A^x + tJ, V = const, (20) 
A< = const, (21) 
w = Mw + V. (22) 

Taking into account that TT~^ is antisymmetric and THH^^T^^ is sym- 
metric part of M. ([T9|) , correspondingly, we get from (21) 

t(t)T"^(t) = const, (23) 
T{t)H{t)H-\t)T-\t) = const. (24) 

Relation (^) yields the system of three ordinary differential equations 
for the functions a(t), /5(t), 7(t) 

a + P cos 7 = Ci, 

/3cosasin7 — 7sinQ; = C2, (25) 
$ sin a sin 7 + 7 cos a = C3, 

where Ci,C2,C3 are arbitrary real constants. Integrating the above system 
we obtain the following form of the matrix T{t): 

Tit) = C,fC2, (26) 

where Ci,C2 are arbitrary constant 3x3 orthogonal matrices and 

- cos s cos bt sin s cos s sin bt 
T=\ sinbt cosbt | (27) 

sin s cos bt cos s — sin s sin bt 
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with arbitrary constants b and s. 

The substitution of equahty (|26|) into with subsequent differentiation 
of the obtained equation with respect to t yields 

C^^f-^fC2L + L + LC^\f-^)fC2 = 0, (28) 

where L = HH~^, i.e. k = hi/ hi, i = 1,2, 3. From (pSD we have 

li = const, i = 1,2, 3; 
b {h — h) cos a2 sin 72 = 0, 

b {h — h) (— sin 02 sin /32 + cos 02 cos /92 cos 72) = 0, (29) 
b {h — ^3) (sin ^2 cos P2 + cos 0:2 sin P2 cos 72) = 0, 

where a2,/32,72 are the Euler angles for the orthogonal matrix €2- Thus we 
obtain the following forms of hf. 

hi = Ciexp{lit), Ci = const, U = const, i = 1,2,3. (30) 

From (^) we get the possible forms of b, li and €2'- 

(i) 6 = 0, h,l2,h are arbitrary constants, 

C2 is an arbitrary constant orthogonal matrix; 

(tt) MO, ll = l2 = h, 

C2 is an arbitrary constant orthogonal matrix; (31) 
(m) by^O, h = l2 7^ h, 

(El cos 9 — ^isin^^ 
-£162 

62 sin 9 82 cos 9 

where 61,62 = ±1, and 9 is arbitrary constant. We do not adduce cases 
6 7^ 0, h ^ I2 = h and 6 7^ 0, h 7^ h = h because they are equivalent to 
case (Hi). 

Now the last equation from the system (0)-(|l3) takes the form 

E^o(-.)^|^ + To(t) + E/. = 0. 
i=i '-'•^j '-'•^j i=i 
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Splitting this relation with respect to independent variables uii,uj2,i03,t for 
each class of functions z = z{uj) given in (|l^) yields the explicit forms of the 
functions Foi(a;i), -Fo2(ci;2), F'q^^uj^) and T'o(t) up to the choice of Aj, i = 1, 2, 3 
in (I) 

3 

F,o = 0, To = -Y.h. (32) 

i=l 

We summarize the above-obtained results in the form of the following 
assertion. 

Theorem 1 The Fokker-Planck equation admits separation of variables 
if the drift coefficients B[x) are linear and given by formulae ^^), where the 
matrix M is defined by formulae ([7^, and (|^. 

The coordinate systems allowing for variable separation in the corresponding 
FPE are given implicitly by formulae (|TT]), (^) and (^^, where T[t) is given 



in (|26|), (p^) and (pT]), functions hi{t), i = 1,2,3 are given in (pO]) and 
functions Wi{t), i = 1,2,3 are solutions of system of ordinary differential 
equations (|22|). Further details on explicit forms of the drift coefficients and 
the coordinate systems are given in Section III. 



III. Exact solutions 

Remarkably, for the equation under study it is possible to give a complete 
account of solutions with separated variables. They have the form and 
the separation equations for the functions ip^, (// = 0,1,2,3) read as (||), 
where the coefficients Fai, a,i = 1,2,3 are the entries of the corresponding 
Stackel matrices (pTSl), functions T^, a = 1,2,3 are hsted in (|16D and the 
functions Tq, Fao, a = 1, 2, 3 given in (152). 



The separation equation for the function ipo{t) is easily integrated. The 
separation equations for the functions ipi^uji), {i = 1, 2, 3) are similar to those 
arising from separation of variables in the Helmholtz equation (A3 + co'^)\l' = 
0. The solutions of these equations are well known (see, |T5| and the 
references therein). Below we adduce solutions of FPE (|1]) for each class of 
functions z = z{uj) given in (p!3D. 
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1. Cartesian coordinates 
r 3 



and Ai = — a^, A2 = — A3 = —7^. 
2. Cylindrical coordinates 

^t,^)=exp|Ai|^e-^'^* + A3|^e-^'3t 



(2/1 + /3)t X 



where Jn is the Bessel function [jl6|, ^l]' ^^^1 Ai = — a^, A2 = — n^, A3 



3. Parabolic cylindrical coordinates 



-2 -2 



u(t, cJ) = exp I + A3^e 



- (2/1 + h)t X 



A/x-i/2(±o-^^i) ^-v-i/2(±o-t^2) e 



where a = e^'^^^{2aY^'^, D,y is the parabolic cylinder function |T^, |14 
and Ai = — a^, A2 = —2afi, A3 = —7^. 

4. For the case of elliptic cylindrical coordinates we have two types of 
solutions 

-2 -2 



u{t, uj) = exp Ai^e-2'^* + A3^e-2^3t _ (2/^ + ^3)^ 



2Zi 2/3 
Ce„(u;i, q) ce„(u;2, g) e*'^^'^, n = 0, 1, 2, . . . , 



X 



u{t, Lu) = exp \ Ai^e^^'i* + A 



2/r ' ■■'2/3'' 



'^3 ^-2^3^ 



(2/1 + /3)t X 



Se„(t^i, q) se„(t^2, g) e"^'^, n = 1, 2, 3, . . . , 

where the even and odd Mathieu functions, Ce„, Se„ are the 

even and odd modified Mathieu functions [1^, |18| and Ai = — 4ga^, A2 = 
2g + c„, A3 = —7^, and c„ are eigenvalues of the Mathieu functions. 
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5. Spherical coordinates 



c7^ 



u{t,u) = exp - 3/it| x 

cuyV±(„+i/2)(«M)P„±™(tanhcu2)e*"^'", 

where J^, is the Bessel function, is the Legendre function WM and 
Ai = — a^, A2 = —n{n + 1), A3 = — m^. 

6. Prolate spheroidal coordinates 

u{t,u) = exp - 3/it| X 

Psjf l(cothcJi, -a^Ai) Ps|,'"l(tanhcu2, -a^Ai) e^'"'^^ 

where m is integer, n = 0,1,2, ... , —n < m < n, Ps™ is the spheroidal 
wave function and A2 = Al["l, A3 = —m^. 

7. Oblate spheroidal coordinates 

u{t,uj) = exp - 3/it| X 

Ps|r'(-aancui, -a^Ai) Ps|r'(tanhcu2, a^Xi) e''^\ 

where m is integer, n = 0,1,2, ... , —n < m < n, Ps™ is the spheroidal 
wave function and A2 = Al["l, A3 = —m^. 

8. Parabolic coordinates 

u{t,uj) = exp - 3/it| e''""^ x 

e""'^^ exp(±me2'^V2)iFi(-iA2/4a + (m + l)/2,m + l,Tzae2'^i) x 
exp(±me2'^V2) 1^1(^2/40 + (m + l)/2, m + 1, ^iae^'^^). 



where iFi is the confluent hypergeometric function |T^ and Ai 



-a^, A3 = — m^. 
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9. Paraboloidal coordinates 

u{t, uj) = exp 2]~^~^'^* ~ 3/it| gc„(ico'i; 2aa, X2/2a) x 

gc„(co'2; 2aa, X2/2a) gc„(2to'3 + tt/2; 2aa, X2/2a) 

or the same form with gc„ replaced by gs^. Here gc„ and gs„ are the 
even and odd nonpolynomial solutions of the Whittaker-Hill equation 
JT^ and n = 0, 1, 2, . . . , and what is more, Ai = — a^, A3 = fin- 



10. Ellipsoidal coordinates 

U{t, uj) = exp |Al|^e-2W _ g^^^j ^^rn^^^^ 

where m is integer, n = 0, 1, 2, . . . , — n < m < n, el™ is the ellipsoidal 



wave function |18] and Ai = Unm, A2 = Xnm, A3 = fXr. 



11. For the case of conical coordinates we have two types of solutions 



u{t,uj) = exp|Ai|^e 2'i*-3/it| cjf J±(„+i)(a/cui) x 

Ec™(a;2)Ec^(cu3), n = 0,l,2,..., m = 0,l,...,n, 
u{t,uj) = exp l^il^e"^''* - 3/it| cjf 4(„+i)(a/c<;i) x 

Es:r(^2) Es:r(^3), ^ = i, 2, 3, . . . , m = 1, 2, . . . , n, 

where is the Bessel function, EcJ^ and EsJ^ are the even and odd 
Lame functions 113 and Ai = — a^,A2 = — n(n + 1),A3 = 



n y 



where c™ are eigenvalues of the Lame functions. 



In these equations we suppose that 7^ 0, = 1, 2, 3). Given the condition 
li = 0, the expressions exp(— 2Zji:)/2/j should be replaced by —t. 

Finally, we give a list of the drift velocity vectors B{x) providing separa- 
bility of the corresponding FPEs. They have the following form: 

B{x) = Mx + w, 

where v is arbitrary constant vector and M. is constant matrix given by one 
of the following formulae: 
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I. M = TLT-^, where 




h,h, h are constants and T is an arbitrary constant 3x3 orthogonal 
matrix, i.e. is a real symmetric matrix with eigenvalues h,l2,h- 

(a) hihih are all distinct. The corresponding FPE has solution 1 
only from the above list. The new coordinates cui, are given 

implicitly by formula 

x = TH{t)z{u)+w{t), (33) 



where z^Q) is given by formula 1 from (|T3D, w(t) is solution of 
system of ordinary differential equations (E^ and 



/ cie'i* 
H{t) = C2e'2* I (34) 
V C3e'3* 

with arbitrary constants Ci, C2, C3. 

(b) h = h ^ h- The corresponding FPE has solutions 1-4 only from 
the above list. The new coordinates uji,uj2, ujs are given implicitly 
by (^3|) , where z{lj) is given by one of the formulae 1-4 from 
(|13[) and H(t) is given by ( p4D with arbitrary constant ci, 02,03 
satisfying the condition ci = C2 for the partially spht coordinates 
2-4 from (1131). 



(c) h = h = ^3, i-e. M = lil, where I is unit matrix. The correspond- 
ing FPE has all 11 solutions, listed above. The new coordinates 
uJi,uJ2, UJ3 are given implicitly by formula (^), where z{u) is given 
by one of the eleven formulae ([13|) and H{t) is given by (|3^) with 
arbitrary constants Ci, 02,03 satisfying the condition ci = C2 for 



the partially split coordinates 2-4 from (|l^) and the condition 



Ci = C2 = C3 for the non-split coordinates 5-11 from (0). 
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COS s 
M = 6 Ci I - cos s sin s I C^^ + IJ, 
— sin s 

where / is the unit matrix and Ci is an arbitrary constant 3x3 orthogo- 
nal matrix, 6, s, li are arbitrary constants and 6 7^ 0. The corresponding 
FPE has all 11 solutions, listed above with h = I2 = h- The new co- 
ordinates tui, ti;2, 1^3 are given implicitly by formula (|ll]), where z{uj) is 
given by one of the eleven formulae (|13|), T(t) is given by (p^)-p7|), 
w(t) is solution of system of ordinary differential equations (^) and 

/ ci 
H{t) = exp(/it) C2 
V C3 

with arbitrary constants Ci, C2, C3 satisfying the condition Ci = C2 for the 
partially split coordinates 2-4 from (|13|) and the condition Ci = C2 = C3 
for non-split coordinates 5-11 from (O). 



/ Hh + h + ik - /a) cos 2s) 



M 



Ci 



V 



—b cos s 
i(/3 - h) sin 2s 



6coss 
h 

—b sin s 



|(/3 - Zi) sin 2s 
6 sin s 

5(^1 + ^3 -(^1-^3) cos 2s) 



where Ci is an arbitrary constant 3x3 orthogonal matrix, 6, s,/i,/2 
are arbitrary constants, li 7^ I3 and 6 7^ 0. The corresponding FPE 
has solutions 1-4 only from the above list with h = h ^ h- The new 
coordinates ui,U2,uj3 are given implicitly by formula (pT]), where i*(a;) 
is given by one of the formulae 1-4 from (0), T(t) is given by PB|), (^) 
and (zm) from (|3T|), w(t) is solution of system of ordinary differential 
equations (UT) and 



Hit) 
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C2e'i* 
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with arbitrary constants ci,C2,C3 satisfying the condition ci = C2 for 
the partially split coordinates 2-4 from (13). 



Note that the above obtained solutions can be used as the basis functions 
to expand an arbitrary smooth solution of the equation under study in a 
properly chosen Hilbert space (for more details, see |^). 

The physical analysis of the obtained results seems to be very interesting, 
but the detailed study of this problem goes beyond the scope of the present 
paper. 



IV. R— separation of variables in the Fokker— 
Planck equation 

In this paper, we restrict ourselves to the choice of separation Ansatz in the 
form (^. Generally speaking, the problem of separation of variables includes 
the search of R-separable solutions of the more general form 

u{t, x) = e^(*'^^Vo(t) n (u;a{t, x),X). (35) 

a=l 

In this case we have an analog of the system of equations (^-(^) 

^ , , duJi duJi dR ^ „ „ dR 
J2 F^oiui)-^-^ + +AR + Baj— + 

^ OXj OXj Ot OXa 

dR dR ^ , , dBa 

+ Toit) + Tr^ = 0. (37) 



dXa dXa dXa 

Equations (0)-(H) are not changed. In a way analogous to that used above 
we get from (|36|) the form of the drift coefficients B{x) 

B{x) = M{t){x-w) +'iv -2VR, (38) 
where Ai{t) is given by formula (|19[). 
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The compatibility conditions of the above system of PDEs (|38| ) yield 
Bix2 - B2X1 = -2{a + (3 cos 7), 

Bix3 — B?yj.^ = — 2(/5cosQ;sin7 — 7sina), (39) 
B2XJ, — B3X2 = ■~2(/9 sin a sin 7 + 7 cos a). 

As the functions 52,-63 are independent of t, it follows from these con- 
ditions that roti? = const and the functions Q;(t), /5(t), 7(t) obey the system 
of ODE (gl). Thus the matrix T{t) have the form @. 
Consequently the following assertion holds true. 

Theorem 2 For the Fokker-Planck equation to be R-separable it is nec- 
essary that the rotor of the drift velocity vector B{x) is constant. 

Concluding Remarks 

— * 

It follows from Theorem 1 that the choice of the drift coefficients B{x) al- 
lowing for variable separation in the corresponding FPE is very restricted. 
Namely, they should be linear in the spatial variables xi,X2,x^ in order to 
provide separability of FPE (Q) into three second-order ordinary differential 
equations. However, if we allow for separation equations to be of lower or- 
der, then additional possibilities for variable separation in FPE arise. As an 
example, we give the drift coefficients 

Bi{x) = 0, B2{x) = 0, 53(x) = fia (T^fT^) , 

where B^ is arbitrary smooth function. FPE (0) with these drift coefficients 

separates in the cylindrical coordinate system t.uji = In ^y^xf + x^^ ,uj2 = 

arctan(xi/x2), ^3 = x^ into two first-order and one second-order ordinary 
differential equations. 

For the one-dimensional FPE the choice of the drift coefficients B[x) 
allowing for variable separation is essentially wider [^]. 
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